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The three-dimensional classical dimer model with interactions shows an unexpected continuous 
phase transition between an ordered dimer crystal and a Coulomb liquid. A detailed analysis of the 
critical dimer and monomer correlation functions point to a subtle interplay between the fluctuations 
of the crystal order parameter and the "magnetic" degrees of freedom present in the Coulomb phase. 
The distribution probability of the crystal order parameter suggests an emerging continuous 0(3) 
symmetry at the critical point. 
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Introduction — The Berezinskii-Kosterlitz-Thouless 
transition is an example of critical phenomenon which 
cannot be understood from the conventional Landau- 
Ginzburg- Wilson (LGW) point of view, based on spon- 
taneously broken symmetries and critical fluctuations of 
the associated order parameter Many quantum sys- 
tems in one dimension (ID) and classical systems in 
2D exhibit such type of continuous phase transition, 
which is understood as a proliferation of topological de- 
fects. The recent suggestion |^ that some phase transi- 
tions in D > 2 could also be outside the LGW frame- 
work triggered an important activity, both on the field- 
theoretical Q and numerical ^, ^ sides. 

The classical dimer model on the cubic lattice ^ 
is one of the simplest 3D models with a transition which 
does not seem to fit in the LGW picture ^ , although the 
field theory of the critical point has not yet been eluci- 
dated. Configurations of the model are hard-core dimer 
coverings of a cubic lattice |^ . Each dimer occupies two 
sites and each site belongs exclusively to one dimer. In 
addition, each configuration has an energy equal to -1 
times the number of square plaquettes with two parallel 
dimers [||. At low temperature T, this interaction favors 
the six-fold degenerate regular (crystalline) arrangements 
of dimers. At high-T, the lattice symmetries are restored 
(dimer liquid), but dimer-dimer correlations decay alge- 
braically with distance (with the same r~^ dependence as 
the interaction energy of two magnetic dipoles) Q . This 
so-called Coulomb phase has also been found in other 
3D models and is relevant to some pyrochlore com- 
pounds (see for instance Rcf. pd||). 

The most intriguing properties of this dimer model 
concern the phase transition between the low-T crystal 
and the Coulomb phase. Previous simulations ^ indi- 
cated that the transition is continuous (second order). 
This is already in contradiction with the LGW approach, 
which generically predicts a first order transition in the 
present case. The main question, still open at present, is 
to understand what are the fields needed for a descrip- 
tion of the long-distance properties at the critical point. 



and what is the associated field theory. In this Rapid 
Communication, we report on numerical investigations of 
correlation functions and order parameters in the vicinity 
of the critical temperature T^. We show that the dipolar 
character of the dimer correlations is gradually lost when 
approaching the transition, and find an unexpected prob- 
ability distribution of the order parameter, with a possi- 
ble emerging continuous symmetry. Related to rotations 
of a 3-component order parameter, this 0(3) symmetry is 
not present at the microscopic level. Terms which make 
the symmetry group discrete at the lattice scale (a dimer 
can only take six orientations) may be irrelevant in the 
renormalization group (RG) sense at the critical point, 
but clearly relevant in the ordered phase. This finding 
is a step toward a description of the transition in the 
continuum limit, and establishes a tight connection with 
previous models (quantum and classical ^) discussed 
in the context of non-LGW transitions. 

Dimer correlations: from dipolar to spin-spin form — 
In Coulomb phases, it is useful to introduce a "magnetic 
field" Q defined on each bond by Bij = —Bji = 1 if a 
dimer sits between sites i E A and j E B {A and B are 
the two sublattices), and Bij ~ otherwise. When con- 
sidered as a vector field B, its lattice divergence satisfies 
V-B = {—y = ±1, as each site is occupied by exactly 
one dimer. Coarsegraining over a large number of sites, 
the contributions of both sublattices cancel out to give 
V • B = 0. In the continuum limit, B is thus written 
using a vector potential: B = V x A. At high-T and 
long distances, this zero-divergence constraint captures 
the essential aspects of the dimer hardcoreness 1^. The 
Coulomb phase is then described by the simplest action 
compatible with the gauge invariance A — > A + VA (A is 
an arbitrary function), that is = / rf^r(V x A)^(r). 
Here the rigidity K^^ depends on T. This effective action 
leads to dipolar long-distance correlations H : 
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K ^ can be measured numerically through LK 
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(((/)^)^), the fluctuations of the total flux 0^ = 
/ dxdyB^{r = [x,y,zo]) going through the sample (in- 
dependent of zq). Ref. ^ showed that K^^ vanishes con- 
tinuously when approaching the transition at T^. ~ 1.675, 
as magnetic field fiuctuations are costly when approach- 
ing the crystal. Although the dipolar form dominates 
at high-T, we find that this contribution to the dimer- 
dimer correlations (Eq. [l|) progressively disappears when 
approaching Tc- In Fig.^A, the correlations {df^dlY for 
r = [i, j, 0] are multiplied by |rp to highlight the dipolar 
contribution, and its strong reduction when going from 
T = 4 to T = 1.68 (just above Tc). This is expected since 
K^^ vanishes when approaching Tc- However, the "mag- 
netic" action does not account for an important qualita- 
tive feature. Setting a — (3 — x in Eq. |l|, the prefactor a 
of the dipolar contribution a/r^ , depends on the lattice 
direction (for instance a — Q for r = and a = 1/2 

for r = [z,i,0]), whereas the actual correlations close to 
Tc appear to be almost isotropic (see Fig. ^[A). 

Eq. |l| misses an important physical ingredient: the 
crystal-like dimer correlations associated to the order pa- 
rameter of the symmetry breaking phase. At each site, we 
define a three-component vector n by n"(r) = ±(—1)''° 
if the dimer sitting in r points in the direction ±a, and 
otherwise. This insures that n is "ferromagnetically" 
correlated in the crystal phase ((n) = ±e^,±e^,±e^ 
in the six ground-states). Making the assumption that 
long-distance fluctuations of n are described by a three- 
component "spin" model [i.e. absence of coupling be- 
tween real space r" and internal spin indices, as in 
0{n) models), the angular dependence of the correlations 
become (n"(0)n''(r))'^ ~ 5°"^ . It turns out that in a wide 
temperature range around Tc, the connected dimer-dimer 
correlations {d°' {Q)d^ {r))Y are well approximated by a 
dipolar contribution, plus a "spin-like" contribution: 
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Taking a = x, (3 = y and r = [i,i,0] (Fig. |i|A), the 
"spin" term vanishes and only the dipolar part remains, 
whereas for a = P = x and r = [i,i,i] the dipolar part 
vanishes and one is left with g{r). The latter is displayed 
in Fig. ^3, where strong and slowly decaying correlations 
are found close to Tc. At Tc, g{r) is very likely algebraic 
and compatible with ^ l/r, although it is difficult to give 
an error bar on the exponent. 

A more quantitative analysis can be done on the 
monomer correlator Z{r,T), defined as the ratio of the 
partition function with two fixed test monomers at dis- 
tance r, by the bare partition function. In the Coulomb 
phase, Z{r,T) converges to a finite value when r — > cx) 
@. As shown in Fig. |, Z(r, T)/Z(oo, T) is weU fitted by 
a function F{x) of the single parameter x — r\/T — Tc, 
suggesting the existence of a length ^(T) ~ (T — Tc)^2 , 
diverging at Tc with an exponent v compatible with 
0.5. By construction F(x oo) = 1, and we ob- 
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FIG. 1: (color online) A: Dimer correlations {dgd"^) as a func- 
tion of |r| and multiplied by |r|"^ for a few selected T above 
Tc (system size L = 64). For this orientation, the spin contri- 
bution vanishes. B: Current correlations (Jo Jr ) of the ICLM 
as a function of |r| and multiplied by |r|'^, at T = 0.5 (in 
the Coulomb phase) and Tc = 0.33305 (system size L = 64). 
Correlators at Tc have been multiplied by 50. Simulations of 
the dimer model (ICLM) were performed with the algorithm 
of Ref. (Ref. 0). 
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FIG. 2: (color online) Dimer correlations {dod^)". A: T ~ Tc 
and just above, for different orientations of r (L = 128). B: 
A few selected T above Tc {L — 64). The lattice direction 
r = [i, i, i] is such that the dipolar contribution vanishes. 



serve F{x) ^ 1/x for x ^ 1, indicating that for dis- 
tances r ^ ^(T), Z(r, T) decays in a way consistent with 
~ l/r. At distances large compared to the lattice spac- 
ing, but small compared to some (diverging) correlation 
length ^(T), the correlations are thus qualitatively not 
of Coulomb type, although T > Tc- Below ^, the cor- 
relations (of both dimcrs and monomers) are compatible 
with an exponent 77 close to zero. Noteworthy, these val- 
ues are those obtained from independent thermodynamic 
measurements 

It is instructive to compare these results with the cor- 
relations of a simpler model: the "integer current loop 
model" (ICLM) |l|, |l|l. There, the allowed configu- 
rations are integer- valued currents Jy- = —Jji € Z de- 
fined on each bond (ij) of the cubic lattice (analogous 
to the dimer magnetic field B), with the constraint that 
the net current arriving at each site is zero: divJ — 0. 
The energy of each configuration is E — ^ S{i j) (Jij)'^- 
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FIG. 3: (color online) Monomer correlation as a function of 
r/^ in log-log scale for a few T just above Tc, assuming ^ — 
(T -ny^ and Tc ~ 1.6743 (L = 128). 



This model also has a Coulomb phase at high-T, with 
large current fluctuations and dipolar (J"J^) correla- 
tions. At low-T, small current loops dominate and even- 
tually Jij = everywhere at T = (a state which does 
not have a dimer analog). Through a standard dual- 
ity transformation ||l5|, the ICLM maps exactly onto an 
0(2) spin model with Villain interaction jl6|. In contrast 
to the dimer model, the phase transition of the ICLM 
(yKLM ^ 0.33305 0) is known to belong to the 3D- 
XY universality class. The Coulomb phase of the ICLM 
corresponds to the ordered phase in the angular variables 
6r, and the dipolar correlations are a direct consequence 
of the Goldstone mode for these dual spins. At high-T, 
this dictionary works also well for the dimer model: 
is the stiffness of a dual 0{2) ferromagnet. The monomer 
correlation Z{r) corresponds to the spin-spin correlation 
(exp (i^o ^ iOr)), and Z(oo) (finite in the Coulomb phase 
and vanishing at Tc) to the square of the order parame- 
ter (magnetization) |(e*^°)p. The crucial difference with 
the dimer model is the absence of spontaneously broken 
symmetry in the low-T phase of the ICLM, and thus the 
absence of an order parameter similar to n. This does 
not only make the critical exponents of the dimer model 
dilferent from that of 3D XY ||], but marks a qualita- 
tive difference concerning the correlation in the Coulomb 
phase close to Tc (for r ^ ^{T)). The current correla- 
tions of the ICLM remain dipolar down to Tc- Both in 
the Coulomb phase and at Tc, they decay as ~ 1/r'^ with 
a geometric prefactor which does depend on the direc- 
tion of r (see the large- |r| behavior in Fig. |^B). The ICLM 
obeys Eq. |l| without the need for an additional "spin like" 
contribution. Only the overall amplitude is reduced close 
to Tc, with l/K{Tc) ^ 1/L (see renormalization factor in 
Fig. |i|A) . The fact that dimer correlations are dominated 
by non-dipolar contributions close to Tc gives additional 
support for this transition not being of 3D XY nature, 
as naively expected from the high-T side. 

Symmetry in the crystal order parameter at Tc — 
The results above point to the importance of fluctua- 
tions of the crystal order parameter n at the transition. 



We therefore looked at the distribution probability of 
N = n(r) close to Tc. For large systems, the 

probability distribution of {N'^ , , ~ 0) is almost 
circular (see the distribution for L = 128 and T = 1.673 
in Fig. I^B). Some small square-hke anisotropy is still de- 
tectable but we argue below that it could be a finite size 
effect and the asymptotic form of the distribution spheri- 
cally symmetric in the thermodynamic limit. To quantify 
this, we define 



C4 = i(cos(4arg[A^- 



iN^]) + cos(4arg[iV^ + iiV^]) 
cos(4arg[iV^ -f iTV^)), (3) 



which measures the 0(3) character of the angular distri- 
bution of iV. If iV is distributed isotropically, 04 vanishes 
in the thermodynamic limit. If N preferentially points 
along the axis directions, one expects C4 > 0. In the 
crystal phase 04 ^ 1 when T — > 00, signaling that the 
direction of the order parameter N is locked along one 
crystal axis (see the 04(T = 1.66) and the associated 
distributions in Fig. |4|A and B). In the Coulomb phase, 
N is typically very small {{N'^) ^ and distributed 

in a Gaussian way, thus with spherical symmetry. The 
non-trivial result is the behavior of C4 at Tc- Fig. ^ 
suggests that C'4(T = Tc, L) ^ when L ^ 00, and thus 
indicate an 0(3) symmetric distribution of the order pa- 
rameter jlTt . 

For T slightly below Tc, Fig. |A reveals that C4(T,L) 
has a minimum for a system size Lq(T). This minimum 
shifts to larger size when approaching Tc from below. 
Lo(T < Tc) thus defines a lengthscale, diverging at Tc, up 
to which the discrete symmetry breaking becomes weaker 
at long distances. 

This is similar to the 3D XY model with Zq>4 
anisotropy ||l^, |T^, where the g— fold anisotropy is dan- 
gerously irrelevant and makes the (2-component) order 
parameter apparently U{1) symmetric in the ordered 
phase, up to some lengthscale diverging at Tc. A crucial 
dilference here is that the order parameter n has 3 com- 
ponents, and a RG analysis of the 0(3) model predicts 
that cubic anisotropics, such as (n^)^ -I- (n^)'' -I- (n^)*, 
should be relevant |l^. The latter anisotropic 0(3) 
model is however qualitatively different from the model 
studied here, as its high-T phase is a featureless param- 
agnet. From this point of view, the classical spin model 
of Ref . H (where a constraint equivalent to divB = is en- 
forced on spin configurations) may be closer to our study 
since it also possesses a high-T Coulomb phase. Dis- 
crete anisotropics seeming irrelevant for the dimer critical 
point, both models would be predicted from the symme- 
try point of view to share the same universality class. 
However, based on the values of the critical exponents, 
this is not the case. 

Discussion — The numerical observations of this 
Rapid Communication emphasize the unconventional na- 
ture of the Coulomb-crystal phase transition. On one 
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versality class of this transition will also become experi- 
mentally relevant if a Coulomb-crystal phase transitions 
was observed in nature, as in a frustrated (ice-like) mag- 
net for instance. 
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FIG. 4: (color online) A: Spherical character of the order pa- 
rameter N distribution (C4) as a function of l/L for different 
T. The data for T = 1.6 have been divided by 2. B: Probabil- 
ity dibUibuLiun r'(iV'^,iV", l A^^ I < 0.03) as a fimLllun uf A^'^ 
and i N^ for Bovoral tomporaturo and eyetom eiaoE. The range 
is A^ r^' £ [-0.3, 0.3] fui L < 128 and [-0.2, 0.2] fui L - 128. 



Electronic address: 
V. L. Berezinskii, Sov. Phys. JETP 32, 493 (1971); J.M. 
Kosterlitz and D.J. Thouless, J. Phys. C 6, 1181 (1973); 
J. M. Kosterlitz, J. Phys. C. 7, 1046 (1974) 
See e.g. L.D. Landau and E.M. Lifshitz, Statistical 
Phystcs 3rd ed., (Pergamon, New York, 1980) 



[2: 

[3] T. Senthil et al, S cience 303, 1490 (2004); |Phys. Rev. B 



70, 144407 (2004) 



73. 134402 (20061 



71. 125102 (2 005 



|4| D.L. Bergman, G.A . Fiete and L. Balents, Phys. Rev. B 
O.I. Motrunich and T. Senthil, ibid 
P. Ghaemi and T. Senthil, ibid 



me presence 01 a aivergmg lengmscaie dciow wnicn 



side 
corr 

by the sole high-T "magnetic" point of view. From the 

nrdprpd sidp wp pvbibitpd a Inngt.hsralp bplnw whirh thp 



054415 (2006)1; K. Gregor and O.I. Mo trunich, ibid 76 



ilations are not (Jouiomb-iikc cannot be interpreted 



174404 (2007)^ T. j Srover and T. Senthil, |Phys. Rev. Lett 
98, 24720272007)1; O. I . Motrunich and A. Vishwanath 



preprint arXiv:0805.1494 



[51 O.I. Motrunic h and A. Vishwanath, Phys. Rev. B 70 



rryskal nrdpr paramptpr sppms t,n fliirtnatps in an iinPY- 

pected 0(3) symmetric way. Besides the difficulty to 
build a LGW theory from the corresponding order pa- 
3ters tnose are clear qualitative mdications tliat 



075104 (2004)1 

|6| A.W. Sandvik, IPhys. Rev. Lett. 98, 227202 (2007)1 ; R.G 

Melko and R.K. Kaul, ibid [LOO, 017203 (2008l|^ 
[71 A. Kuklov, N. Prokof'e v and B. Svistunov, phys. Rev 



individual (magnetic or crystallmej descriptions cannot 
capture the physics involved at the transition. This 
points to the existence of a critical regime where the 
interaction between the "magnetic" degrees of freedom 
at play in the Coulomb phase and the "spin" degrees of 
freedom which order in the crystal dominates the long 
distance physics. This phenomenology is similar to that 
of a non-compact CP^ (NCCPl) model ||, |^, which pre- 
cisely couples spin and gauge degrees of freedom. The 
NCCPl theory is characteriypd b}^ two cnntinunus sym- 
metries" STJ{?) for spin rntatinns, and U{^) fissnciated 



"181 

[9] 
[101 



Lett. 93, 2304 02 (2004) 



D. Huse et a/. |Phys. Rev. Lett. 91, 167004 (2003) 



Phys. Rev. Lett. 97, 030403 (2006) 
d State Gommun. 10. 967 (1972 



F. Alet et al. 
J. Villain, So: 
Isakov et al., 
Banerjee et al. 
T.S. Pickles, T. E. Saunders and J.T. Ghalker, preprint 



Phvs. Rev. Lett. 



93. 



167204 (2004) 



S.V. 
A. 



Phys. Rev. Lett. 100, 047208 (2008) 



[11] 



arXiv:0708.3791 



et al. 



Phys. Rev. Lett. 95, 217201 (2005 



S. BramweU et al., [Phys. Rev. Lett. 87, 047205 (2001)1 ; 
G.L. Henlev, |PIw"s. Rev. B 71, 014424 (2005)| ; S. V. Isakov 



[121 A.W. S andvik and R. Moessner, [Phys. Rev. B 73, 144504 



(2006) 



to the magnetic flux conservation. In our case, the flux 



con^ jCTvation is implemented at the lattice level, and 



we 



[131 M. Wallin et al, [Phys. Rev. B 49, 12115 (1994) 
[141 F. Ale t and E.S. S0rensen, jPhys. Rev. E 6 7, 015701 



(2003) 



argued that continuous "spin" rotations emerge at the [15I J. V. Jose et al, ^hys. Rev. B 16, 1217 (1977) 



critical point. This is similar to the 2D quantum spin 
system studied in Ref. ||, but there the spin rotation sym- 
metry is explicit and a gauge flux conservation is argued 
to e merge dynamically at the critical point. Doopito all 



[161 J- Villain, J. Phys. (Paris) 36, 581 (1975) 
[171 Since T — 1.675 may be slightly below the actual Tc, 

Ci{T = 1.675) may start to increase for L > 128 
[181 J. Lou, A. W. San dvik and L. Balents, [Phys. Rev. Lett 



thes e common ingrodionto between NCCPl thoorioD and 



99, 207203 (2007) 



our idimer model, there is a clear mismatch between our 



[191 J. M. Garmona A. Pelissetto. and E. Vicari, Phys. Rev 



critical expon(ints 'and th6s6 iittributfid so f;i,r tO NCCPl [20] 
critical points ||, ^. So, the actual field theory describ- 
ing the dimer model transition remains to be uncovered. 
Besides its theoretical importance, the issue of the uni- 



B 61, 15136 (2000) 



F. Albuquerque et al., J. Magn. Magn. Mater. 310, 
1187 (2007); M. Troyer, B. Ammon and E. Heeb, 
Lect. Notes Gomput. Sci. 1505, 191 (1998). See 



http : / / alps . comp-phys . org 



